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Abstract 

In this work, a spin I relativistic particle described by a generalized potential containing both 
the Coulomb potential and a Lorentz scalar potential in Dirac equation is investigated in terms 
of the generalized ladder operators from supersymmetry in quantum mechanics. This formalism 
is applied for the generalized Dirac-Coulomb problem which is an exactly solvable potential in 
relativistic quantum mechanics. We obtain the energy eigenvalues and calculate explicitly the 
energy eigenfunctions for the ground state and first excited state. 
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I. INTRODUCTION 



Using first the separation in polar coordinates system for the generalized Dirac equa- 
tion, in three dimensions, the supersymmetry (SUSY) in quantum mechanics (QM) is inves- 
tigated. We consider a spin | particle in the Coulomb potential, K = — > an d a Lorentz 
scalar potential, added to the mass term in the Dirac equation which may be interpreted 
as an effective position dependent mass. If the scalar potential is assumed to be created by 
the exchange of massless scalar mesons, it has the form V s = — ^+ 

The (3+1) and (1+1) dimensional Dirac equations with both scalar-like and vector-like 
potentials have been well known in the literature for a long time [l]. Exact solutions for 
the bound states in this mixed potential can be obtained by the method of separation of 
variables 2|-|5| and also by the use of the dynamical algebra SO (2, 1) [6(. In a recent paper 
the solution of the scattering problem for this potential has been obtained by an analytic 
method and also by an algebraic method by Vaidya and Silva Souza 

We do not think that the solution of the relativistic Coulomb potential with position- 
dependent mass in the Coulomb field has been correctly solved. 

Recently exact solutions have been found for fermions in the presence of a classical back- 
ground which is a mixing of the time-dependent of a gauge potential and a scalar potential 
[9]. Also, exactly solvable Eckart scalar and vector potentials in the Dirac equation have 
been investigated via SUSY QM (lOj, the 5- wave Dirac equation has been solved exactly 
for a single particle with spin and pseudospin symmetry moving in a central Woods-Saxon 
potential JjJ, and it was shown that some other cases for which the Dirac equation with 
classical potentials of vector and scalar natures with spherical asymmetry can be solved 



exactly 



12|. 



In this work we consider an alternative calculation for the energy spectrum and eigenfunc- 
tions of the Dirac equation via a Schrodinger-like wave equation, for the vector and scalar 
pot entials recently studied by Alhaidari [8]. The connection between the Johnson-Lippmann 



ij, [jjl and the generalized Johnson-Lippman operators 
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16| in a Dirac-Coulomb problem 



via SUSY QM has been studied in 3-dimensional space. In 13J relations between various al- 

n 

gebraic approaches via radial equation are pointed out, and in [14j the bound eigenf unctions 
and spectra of a Dirac hydrogen atom have been found via su(l, 1) Lie algebra. 

Our calculation generalizes the results previously obtained for the relativistic Coulomb 
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problem. Instead of a direct generalization of Sukumar's calculation we use the equivalent 



approach developed by Fukui and Aizawa [19| and also by Balantekin [20|, |2l| where shape 



mvanance 



22| plays an important role. The Coulomb problem in non-relativistic quantum 
mechanics has been obtained by the use of SUSY shape invariances j^]. The special case 
of the Dirac-Coulomb potential has been treated recently via SUSY quantum mechanics by 



one of the authors 



28|. 



This work is organized as follows. In section II, from the time independent Dirac equation 
for a potential in terms of the Coulomb potential and a Lorentz scalar potential we construct 
the radial equations. In section III, we have achieved the diagonalization of the matrix that 
appears in the interaction term of the radial Dirac equation, supercharges, supersymmetric 
Hamiltonian and the ground state energy eigenvalue and eigenfunction. In section IV, using 
the shape invariance properties we deduce generalised ladder operators to build up the bound 
state energy eigenvalues and eigenfunctions. The conclusion is presented in section V. 

II. TIME INDEPENDENT DIRAC EQUATION 

The time independent Dirac equation with vector V(r) and scalar potentials Vs(r) may 
be written in the form 

HV = Eif), (1) 

where the Hamiltonian is given by 

H = Pl ®a-p + p 3 ® l 2 x 2 (M + V s (r)) + V(r)l 4X 4, (2) 

and we have used a direct product notation in which pi and Cj, {i = 1,2,3) are the Pauli 
spin matrices obeying [pi, erj]_ = 0, with % — 1 = c. 

The Hamiltonian (j2J) commutes with the total angular momentum 

J = L+ l -a. (3) 
It also commutes with the Dirac operator 

K = p 3 (l + a-L), (4) 
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and with the inversion operator 



P = P*I, 



(5) 



where / inverts the spatial coordinates and the momenta. 

The symmetry group of this Hamiltonian (j2J) is characterized by the two vector invariants 
J e K, precisely as in the nonrelativistic case. A complete set of mutually commuting 
operators is H, J 2 , J 3 , K, P which have the simultaneous eigenvectors | E,j, m, k > with 



■J 2 \E 1 j 1 m,K> = j(j + 1) | E, j,m,K > 
J3 I E, j,m, k > = m I E, j, m,K> 
K I E, j, m, k > = — k I E, j, m, k > 



(6) 



where j half integral, m = — j, ■ ■ • j and k — ± y j + ~ 

To separate variables we introduce two component spinors which are eigenf unctions of 
J 2 , J Z ,L 2 , S 2 and are of two types 



for j = t + i, k = j + \ > and 



21+1 



Yi 



\ 



£,m- 



\ { 21+1 



(7) 



3,m 



21+1 



V 



for j = l-l,k = -(j + ±). 

In the above basis one verify that 



Y, 



t,m- 



1 e,m+ l 



21+1 



(8) 



(9) 



a ■ ni 



(10) 



where n 



and 
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(l + a-L)cf>% = -k<t> < g, (11) 

where k — ± (j + |) for j — 
Next, we put 



where 



/ ±1 

r ^jm 



(12) 



<^m = 0jm\ (13) 

for j = £ ± | . The operator a ■ n when expressed in this representation, as was done in (FTU]) . 
has the property of reversing the sign of k. The result which can be easily established by 
means of this operator relation 

[1 + <?• L,a- n} + = 0. (14) 

Next, using this relation we have 

Km = -km. (15) 

Further, using the relations 



cr . p f = + _ \ a . n( p ( 16 ) 

7^ If \ CLf f 



and 



- - -/( r ) ^ * f d f k f\^e m 



we get the radial equations 



^+*G ij -(E + M + V s (r) + V(r))F ij = 0, 



dFJy k 



dr r 

where k is the eigenvalue of the Dirac operator K 



-F £j + (E-M -V s (r) + V(r))G ej =0, (18) 



5 



Defining 



A = ±Kp 3 - rV s {r)pi - irV{r)p 2 , 
k- p = Mpi +iEp 2 , 

the Dirac radial equation for the hydrogen atom may be written in the form 



(19) 



d a r _; 

— + --k-p 

ar r 



$ = 0. 



(20) 



Next we show that as the Dirac equation becomes a Schrodinger-like wave equation, with 
the vector and scalar potentials. 

III. ENERGY SPECTRUM AND EIGENFUNCTIONS VIA SUSY 

Now, let S be the operator which diagonalises the matrix that appears in the interaction 
term. The matrix A may be diagonalized by using the result 



where 



and 



S-'AS = Ap 3 , 



X+\k\ 



(21) 



(22) 



(23) 



Defining 



we get 



d A 

T + -P 3 

ar r 



k ■ p 



$ = 0, 



(24) 



(25) 
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where 



2-\{\ K | +A)L = (| k | +A) 2 (M + £7) - (A x - A 2 ) 2 (M - £) 
a 

2-A(| k | +A)A;+ = (| k | +A) 2 (Af - £) - {A 1 + A 2 f{M + E) 



c 



A; 3 = . (26) 



Thus, we get 



a 2 



h '- k+ ~Y 2 



12 i\,r2 



k 2 M 2 =k 2 = M 2 -E 2 (27) 



where 



M 2 a 2 = (EA l + MA 2 ) 2 - fc 2 A 2 , (28) 

which provides a relation between the SUSY and an additional constant of motion in the 
Dirac- Coulomb problem will be reported separately. 
It may be noted that Eq. ff25l) can be rewritten as 

A+G = LP, A~F = -k+G (29) 

where 

A ± = ± d + X_EA 1 + MA 2 
dr r A 

These relations are similar to the relations between the two components of the eigenfunctions 



of a supersymmetric Hamiltonian 17 ] 



n 



[Q,Q t 



V 



= QQ f + Q f Q 

A + A~ \ _ I U- 
A~A + ) \ H + 

A + 




[Q,H]_ = 0=[Q t ,^]. 



(31) 



where Q and are the supercharges satisfying the anticommutation relations of the SUSY 
algebra. The supersymmetric partner Hamiltonians 1-L± satisfy the following eigenvalue 
equations: "H±V± = S± V± • Next, Eq. (1251) may be written as 



2 2 

A'A+G = ^M 2 G, A + A~F = ^—M 2 F (32) 
A 2 A 2 

which shows that every eigenvalue of A + A~ is also an eigenvalue of A~ A + except when 
A~F = 0. This condition gives the lower component for the "ground state" wave function 

F = r x e x (33) 

and 



Eo = ,(, Af + Al A\A\ y A 1 A 2 

M { X 2 + A 2 (A 2 + A\) 2 ) ~X 2 + A 2 ' 1 ' 

where E is the ground state energy eigenvalue, for a 2 = 0. 



IV. EXCITED STATES VIA SHAPE INVARIANCE CONDITION 



The shape-invariant SUSY partner potentials 22|-|25| are similar in shape and differ only 
in the parameters that appear in them. More specifically, if V_(x;ai) is any potential, 
adjusted to have zero ground state energy = 0, its SUSY partner V + (x; a\) must satisfy 
the requirement 



V+(x;ai) = V-(x;a 2 ) + R(a 2 ), a 2 = /(c&i), (35) 

where a\ is a set of parameters, a 2 a function of the parameters a\ and R(a 2 ) is a remainder 
independent of x. Then, starting with V\ = V^(x; a 2 ) and V 2 = V + (x; ai) = V\(x\ a 2 ) +R(a 2 ) 
in ( 135|) . one constructs a hierarchy of Hamiltonians 

H n = + V^x; a n ) + ^: =2 R(a s ), (36) 

where a s = f s (ai), i.e., the function / applied s times. In view of Eqs. ( 1331) and ( 1331) . we 
have 
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H n+ i = + a n+1 ) + S-+ 2 x i2(a s ) (37) 

= + V+(x; a n ) + Z: =2 R(a s ). (38) 

Comparing ( 1361) . (1371) and (}3~8l . we immediately note that H n and i/ n +i are SUSY partner 
Hamiltonians with identical energy spectra except for the ground state level 

= ^ s=2 R{a s ) (39) 

of H n , which follows from Eq. ( 136]) and the normalization that for any V^(x; a) , E_ = 0. 
Thus, we get 

n+l 

E n l =Er l = ... = E^ +l = Y J R{ds), n = l,2,... (40) 

and 



ip[ n) oc A+(a;; ai)A+(x; o 2 ) . . . A+(x; a^ip^^x; a. 



(o) 



ra+l, 



(41) 



Equations ( )40l) and ( 14T|) . succinctly express the SUSY algebraic generalization, for various 



Wn+l 



tpi oc e" 



shape- invariant potentials of physical interest [22|, |24j, of the method of constructing energy 
eigenfunctions (ipj^l) for the usual ID oscillator problem. Indeed, when a± — a 2 — . . . — 
a n = a n+1 , we obtain ip^l oc (a+) n ^ 0) , Af = ■ ■ ■ A+ = a+, tp^ 
where u is the angular frequency. 

The shape invariance has an underlying algebraic structure and may be associated with 
Lie algebra 21] . Now, we present our own application of the shape invariant method outlined 
above. Consider next the supersymmetric partner Hamiltonians given by 



U- 



A + A 
dr 2 



d 2 A(A - 1) 



+ V-(r, A) 



2(EA 1 + MA 2 ) (EA 1 + MA 2 ) 



A 2 



(42) 



and 
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U + = A 'A + = -— + A ( A+1 ) _ 2 (^i + MA 2 ) + (EA 1 + MA 2 f 



dr 2 

■£+V + (r,A). 



A 2 



These SUSY partner potentials are shape invariant, since 



V+(r, A) = VL(r, A + 1) + i?(A + 1), 



where 



i?(A + l) 



+ MA 2 ) 2 {EA l + MA 2 ) 



A 2 



(A + l) s 



Hence 



>T(A)A + (A) = A+(A + 1)A~(A + 1) + i?(A + 1). 



Following the approach of Fukui and Aizawa 19] and also of Balantekin 20|, |21] 
the following ladder operators 



B-(X) = rt(A)A-(A), S + (A) = (B-)t(A), 
where T(A) is a translation operator defined by 



T(A) = e&, 



with T^(A) = e ex , and get 



[5-(A),fi + (A)] = J R(A). 



It is easy to verify that 



n 

B + ) (A)] = £i2(A + 0(B + ) (A). 



i=i 



Hence the energy eigenvalues of the Hamiltonian H_ are given by 
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= ]T/?(a + 



i=l 



(51) 



(A + n) 2 ' A 2 ' " tv ~ / ' E n 

Thus, from equations fl27|) . fl28|) . fl32|) and fl5T|) . the energy eigenvalues associated to the 
lower component F n are given by 



E (n) 



M 2 



\ 1 H 1=^ 



n = 0,1,2, 



Solving the last equation we get 



(52) 



En 

M 



A X A 2 



in 2 + Ax 



± 



A\A\ 



+ 



h 2 -A\ 
n 2 + A x 2 



(53) 



n 2 + V 

where n = n r + A and n r = 0, 1,2, • ■ -. The corresponding eigenf unctions associated to the 
lower component are given by 



F n (r) = (B + ) n F (r) (54) 

where F Q is given in Eq. ( 133]) . 

From the previous equation we can determine the first excited state of the eigenfunction. 
Using 



£+(A) = T(X)A + (X) 
and according to the equations (l30~l) . ( 1471) and ()48l) found 



(55) 



(56) 



Hence, for the first excited state we get the corresponding eigenfunction associated to the 
lower component 



F 1 (r) = B + F Q (r) 



d A EAl+MA 2 \ 8_ x (E A 1 +MA 2 )r 

— -i e»r e * 

dr r X 



(57) 



so, 
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. > „ r (l EA\ + MA 2 \ X-H (EA 1+ MA 2 )r 

F '' r » = iV '(" A(A + 1) J r e ™ ' PS) 
where iVi is a normalization constant for the first excited state. 

V. CONCLUSION 

In the non-relativistic supersymmetric quantum mechanics formalism, ladder operators 
are constructed to obtain the complete set of the bound state energy spectrum and eigen- 
functions for a relativistic potential given by a Coulomb-like and a Lorentz scalar potential 
in the Dirac equation; one adopts polar coordinates. These new generalized ladder operators 
obtained via supersymmetry shape-invariant Hamiltonians in quantum mechanics can be re- 
duced to the ones of the Coulomb potential, which provides us the exact relativistic energy. 
These Hamiltonians, "H_ and "H+, are SUSY-partner Hamiltonians with identical energy 
spectra, except for the ground state level, satisfying two Schrodinger-like wave equations, 
with the vector and scalar potentials. 

The vastly simplified algebraic treatment within the framework of the primary SUSY 

functional space 



26]. The primary 



algebra, in terms of conserved operator (Runge-Lenz operator) in the ful 
for the nonrelativistic Coulomb problem with spin has been investigated 
SUSY with Johnson-Lippmann [29[ operator for the Dirac-Coulomb problem is not possible 

13]. 

Our results on the connection between the generalized Johnson-Lippmann operator and 
SUSY for the bound eigenstates of the generalized relativistic Coulomb problem with posi- 
tion dependent mass will be reported separately. 
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